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We report on a lattice simulation result for four-dimensional jV — 1 SU(2) super Yang-Mills 
theory with the dynamical overlap gluino. We study the spectrum of the overlap Dirac operator at 
three different gluino masses m = 0.2, 0.1 and 0.05 with the Iwasaki action on a 8 3 x 16 lattice. 
We find that the lowest eigenvalue distributions are in good agreement with the prediction from the 
random matrix theory. Moreover the mass dependence of the condensate is almost constant, which 
gives a clean chiral limit. Our results for the gluino condensate in the chiral limit is (\j/y)rr. = 
0.63(12), where ro is the Sommer scale. 
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1. Introduction 

The ,JV = 1 global supersymmetry in four dimension is phenomenologically interesting |jl]]. 
The ,jV = 1 super Yang-Mills theory (SYM) is a simple model containing only gluon and gluino. 
The model is nontrivial since it is expected to have interesting nonperturbative properties like con- 
finement and discrete chiral symmetry breaking. An important physical quantity is the gluino 
condensate which is an order parameter for the spontaneous breaking of the discrete chiral sym- 
metry. If we extend the model to include matters, the Witten index [||] can vanish implying that a 
dynamical supersymmetry breaking can occur with the condensate as an order parameter 

Lattice gauge theory could be a means to study the nonperturbative aspects of supersymmetry. 
However the lattice discretization explicitly breaks supersymmetry and the supersymmetry violat- 
ing operators at the quantum level give rise to the fine tuning problem in general. The jV = 1 SYM 
is exceptional, since the mass term is the only relevant supersymmetry violating operator, which 
can be controlled by taking the chiral limit [§]. The model has been studied on the lattice with the 
Wilson fermion ^ and the domain wall fermion [||] so far. Still these approaches have difficulties 
in obtaining precise chiral limit due to an explicit breaking of chiral symmetry or a large residual 
mass. Even on the lattice, the question whether the R-symmetry of the model is spontaneously 
broken or not, remains unclear. 

Through this work we provide a numerical evidence for nonzero gluino condensate in JY = 
1 SYM by a nonperturbative lattice study. The overlap fermion ^ is a natural choice for the 
model because of its exact lattice chiral symmetry In this study, we use the same strategy as 
our previous studies in QCD [^], where clean chiral extrapolation for the chiral condensate was 
achieved by dynamical simulation of the overlap fermion and comparing the Dirac spectrum with 
the Banks-Casher relation [[R]] and the random matrix theory (RMT) [11]. The key benefit of this 



strategy is that it is free from both of the UV and IR problems when we determine the physical 
value of the chiral condensate, as we describe later. 

2. The model and lattice formulation 

The Euclidean action for ,jV = 1 SYM is 

S = J ^xTr^F^f + hpfD^} . (2.1) 

Here y = ^/"(x) T" denotes the gluino field, which is an adjoint Majorana fermion satisfying y = 
\j/ T C. The action is invariant under the supersymmetry transformation 



and the R-symmetry 



8A^ = -2g\j/Y^e , 8\\f = --a tlv F llv e , (2.2) 



v g-«"a» v : ijf ¥e - ia K , (2.3) 



which coincides with the £7(1) axial symmetry. This U(1)a is anomalous and the axial current 
J ii = YYii 75 V satisfies 

No 2 ~ 

J„ = cS F a ^ v Ff (2 4) 
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The anomaly leaves discrete subgroup Z2n c invariant at quantum level but it can be spontaneously 
broken down to Z2. An order parameter is the gluino condensate (ijfy). No Nambu-Goldstone 
boson appears since Z2n c is a discrete group. 

On the lattice the overlap Dirac operator has no sign problem, since it possesses Ginsparg- 
Wilson relation, 75 hermiticity, and two-fold degeneracy from the adjoint representation. Therefore 
we can write the fermionic action as 

S / = -Ilogdet(Dj v Dov) , (2.5) 

where the overlap Dirac operator D ov in the adjoint representation is constructed from the Wilson 
Dirac operator Dw as 

D m =M (l+D w /^Dl / D w ^j , (2.6) 
(D w f/ X = 5y X 8 ba - K £ [8 y ^V^ + 8 y+fl , x (Vy h Z) T ] • (2.7) 

Vxfi is the adjoint link variable defined as VjSf = 2Tr (Ux^T b U xtl T a ) and the hopping parameter is 
set to K = 8 _j, Mo . We introduce an explicit mass m for gluino in order for a comparison with the 
random matrix models. 



3. Banks-Casher relation and random matrix theory 

The Banks-Casher relation [ |10| ] relates the Dirac spectrum to the chiral condensate. Moreover 
the details of low-lying Dirac eigenvalue distributions are described by the chiral random matrix 



theory (ChRMT) QUO including finite volume corrections to the Banks-Casher formula. 

There are two advantages in extracting the condensate via ChRMT. One is that the Dirac 
spectrum is free from UV power divergences, which we can explicitly see in the extended Banks- 
Casher formula : 

np{X) = lim±(- (W)\m=iX+e + {W)\m=iX-e) • (3-1) 

Another advantage is that the finite volume effects and finite mass effects can be corrected by 
ChRMT. In RMT, the low-lying eigenvalue distributions are characterized by the dimensionless 
parameter /I = mLV. For any given volume V and mass m, the condensate £ is the only free 
parameter that is defined in the infinite volume and massless limits. 

The gluino is an adjoint representation fermion and the corresponding Dirac operator satisfies 
[D,C l K] = 0, where C is the charge conjugation and K is the complex conjugation operator. 
Consequently the universality class for jV = 1 SYM is different from QCD, and it belongs to the 
chiral symplectic ensemble (ChGSE). The partition function is the same as that for chiral unitary 
ensemble, which describes QCD, except that the matrix elements are real quaternion instead of 
complex number. This is equivalently described by partially quenched chiral perturbation theory. 

For later convenience, we present the first eigenvalue distribution of ChGSE for zero topolog- 
ical charge, 

=ce 1 C,}l — — — — — , (3.2) 

£dtnU(2tix) 
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m 


dx x N x 


Mrj 


accept 






r /a 


0.20 


0.25x4 


396 


0.700(29) 


0.977(58) 


0.68747(21) 


4.34(17) 


0.10 


0.20x5 


372 


0.850(15) 


0.971(18) 


0.68752(19) 


5.25(30) 


0.05 


0.25x4 


382 


0.658(26) 


0.996(42) 


0.68795(16) 


5.44(33) 



Table 1: Summary of the simulation, m is the gluino mass, iVtrj is the number of thermalized trajectory, P is 
the plaquette, and ro is the Sommer scale. 



where £ = XLV , jX = rriLV, p, = C, 2 + pi 2 , c is a normalization constant, and X is the eigenvalue. 
The general formula for the higher eigenvalues is given in [12]. We note that the quenched limit 
can be obtained by simply taking the limit \i — > oo. 



4. Results for Dirac spectrum and gluino condensate 

Our dynamical simulation for ,jV = 1 SYM with the overlap fermion is performed on a 8 3 x 
16 lattice with SI/ (2) gauge link fields. For this small lattice, we need to carefully pick up a 
gauge action and the parameters /3 , Mq to avoid the Aoki phase. The /3 is chosen in such a way 
that the lattice has a sufficiently large physical volume, and the Mq is chosen to guarantee the 
locality. The hermitian Wilson-Dirac spectrum for the quenched theory suggests that the region for 
an implementation of the overlap Dirac operator inside two fingers of the Aoki phase is narrower 
in the present case of the adjoint representation than in the case of the fundamental representation. 
Based on a Dirac spectrum survey in the quenched theory, we use the Iwasaki gauge action with 
J3 = 1.05 and M = 1.9. 

The overlap fermion has exact zero modes which are related to the topological charge of the 
gauge fields by the Atiyah-Singer index theorem. The topology plays a significant role on the small 
volume for a comparison to ChRMT. In our simulation, the global topological charge is fixed by 

H 2 

introducing extra Wilson fermions to produce a Boltzman weight det w 2 where Hw = YsD w is 

the hermitian Wilson-Dirac operator with K = & _ l 2M(j ■ This term suppresses the near zero modes of 
Hw and generates a gap, which supports the exponential locality of the overlap Dirac operator and 
reduces the numerical cost for the sign function in the overlap Dirac operator [|3|]. The inclusion of 

H 2 

det i w i should not affect the low energy physics. We select zero topological charge and small 

n w +m tln 

twisted mass of m tm = 0.01. The choice of the twisted mass does not affect the low-lying Dirac 
spectrum in the quenched theory. 

With the Hybrid Monte Carlo algorithm, we generated ensembles with three different gluino 
mass as shown in the table |[ For the low-lying eigenvalue distribution of the overlap Dirac oper- 
ator, 28 lowest eigenvalues are measured at every other trajectory by implicitly restarted Lanczos 
method. Note that the adjoint representation gives two-fold degeneracy, and the 75 hermiticity en- 
forces eigenvalues to form complex conjugate pairs. Therefore there are essentially only 7 different 
eigenvalues among 28. Our main interest is in the lowest eigenvalue, for which the data shows clear 
complex conjugate pairing with two-fold degeneracy. 

The integrated autocorrelation time is estimated by the history of the lowest Dirac eigenvalue, 
which is smaller than 20 trajectories for all three input gluino masses. All statistical errors are 
estimated by the jackknife method with a binsize of 20 trajectories. 
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Figure 2: The lowest eigenvalue distribution is compared with ChGSE. Curves are the predictions from the 
ChGSE for various masses. The condensate determined from (4. 1) for each mass was used as an input. The 
horizontal errorbar (right) reflects the uncertainty of the condensate. 

The relationship between ChRMT and the chiral perturbation theory (ChPT) suggests an iden- 
tification £ = XLV, from which we can extract the condensate. The condensate at each gluino mass 
is determined by comparing the average eigenvalues of lattice data and that of ChGSE. 

poo 

(C)rmt= / dC£ Pl (£,mZV) = (A)iatEy. (4.1) 
Jo 

Figure [j] shows the gluino mass dependence of £. The comparison with ChGSE gives almost the 
same values of the gluino condensate for different input gluino masses. The mass dependence is 
weaker than that in QCD, which can be explained by the absence of Nambu-Goldstone modes, or 
the absence of chiral logarithms. The constant fit gives £ = 0.00393(18). With the Sommer scale 
obtained from the lightest mass (m = 0.05), we obtain (ij/y)^ = 0.63(12) in the chiral limit. 1 

Let us look into more details of the eigenvalue distribution. Figure @ shows the lowest eigen- 
value distribution /?i(£,/i) (left) and the cumulative distribution J () dC,' p\(C,\\i) (right) for each 
mass, where we used the determined condensate value as an input {\i = mLV). The distributions 

'The value l}l^ro = 0.86(5) is slightly larger than that for QCD in our previous study : E 1//3 ro ~ 0.6. 
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Figure 3: The low-lying eigenvalue distribution is shown for each mass. Black curve shows E+ g|r^ 3 , 
where the second term comes from the free theory in 4-dimension. 

show clear mass dependence, and the data at m = 0.1 and 0.2 agree well with the prediction from 
the ChGSE. There is a common tendency that the distribution is narrower than the predicted curves, 
and this tendency becomes clearest for the lightest mass data. Since ChGOE or ChGUE has wider 
distribution compared to ChGSE, ChGSE fits best to the data anyway. This deviation is most likely 
from the higher order corrections of ChPT, which is not considered in RMT. 

Figure 3 shows the low-lying eigenvalue distribution, which can be roughly divided into three 
regions : the lowest edge near X « 0.05, a plateau region between 0.05 < A < 0.25, and an in- 
creasing region for A > 0.25. The overall shape is consistent with the behavior ^tA 3 of the free 
theory plus the condensate 0.00393(18), which is shown as the black curve in the figure. Note that 
the plateau values, which are without leading volume corrections, do not depend on masses very 
much. Since the difference from these values to the obtained condensate value in the chiral limit is 
around 20%, we expect that the systematic error from the higher order effects are less than 20%. 

5. Summary 

We have performed a dynamical simulation of 4d jV = 1 SYM with the overlap fermion for 
the gluino on a 8 3 x 16 lattice. The exact chiral symmetry of the overlap fermion and the predic- 
tion from ChRMT allow us to extract the gluino condensate {xj/yf}^ = 0.63(12) in the chiral limit. 
The lowest Dirac eigenvalue distribution for two heavy masses (m = 0.2,0.1) nicely matches with 
ChGSE. For the lightest gluino mass (m = 0.05) data, the distribution shows a deviation from the 
prediction possibly due to corrections from the next order finite volume effect. The expectation 
value of the lowest Dirac eigenvalue is used to determine the condensate, and the obtained con- 
densate for each mass does not much depend on the mass. This confirms that our strategy works 
very well in the J/ = 1 SYM, since there is no Nambu-Goldstone boson nor chiral log effect. The 
low-lying eigenvalue distribution is consistent with the simple sum of the constant (condensate) 
and the free theory result. 

The only relevant supersymmetry violating operator in the model is the mass term, which 
can be controlled by taking the chiral limit. Therefore it is interesting to see the violation of the 
supersymmetry Ward identity and its dependence on the input gluino mass. 



6 



Lattice study of 4d jY = 1 SYM with dynamical overlap gluino S.-W. Kim 



We thank P. H. Damgaard, J. Giedt, G. Bergner, Y. Kikukawa, and N. Yamada for useful 
discussions. Numerical simulations are carried out on IBM Blue Gene/L at KEK under the support 
of its Large-scale Simulation Program (Nos. 09/10-09, 10-11) and Hitachi SR16000 at YITP This 
work is supported in part by the Grant-in-Aid from the Japanese Ministry of Education, Culture, 
Sports, Science and Technology (No. 20105002, 20105005, 21674002, 23105710), the Grant-in- 
Aid for Scientific Research (No. 20540286, 23244057) from JSPS, and the HPCI Strategic Program 
of Ministry of Education. 

References 

[1] E. Witten, Dynamical Breaking of Supersymmetry , Nucl. Phys. B 188 (1981) 513. 

[2] E. Witten, Constraints on Supersymmetry Breaking, Nucl. Phys. B 202 (1982) 253. 

[3] A. E. Nelson and N. Seiberg, R symmetry breaking versus supersymmetry breaking, Nucl. Phys. B 
416 (1994) 46 [arXiv:hep-ph/9309299]. 

[4] G. Curci and G. Veneziano, Supersymmetry and the Lattice: A Reconciliation?, Nucl. Phys. B 292 
(1987) 555. 

[5] For example, I. Montvay, Supersymmetric gauge theories on the lattice, Nucl. Phys. Proc. Suppl. 53 
(1997) 853 [arXiv:hep-lat/9607035]; R. Kirchner et al. [DESY-Munster Collaboration], Evidence for 
discrete chiral symmetry breaking in N=l supersymmetric Yang-Mills theory, Phys. Lett. B 446 
(1999) 209 [arXiv:hep-lat/98 10062]. 

[6] For example, J. Giedt, R. Brower, S. Catterall, G. T. Fleming and P. Vranas, Lattice super-Yang-Mills 
using domain wall fermions in the chiral limit, Phys. Rev. D 79 (2009) 025015 [arXiv:0810.5746]; 
M. G. Endres, Dynamical simulation ofN=l supersymmetric Yang-Mills theory with domain wall 
fermions, Phys. Rev. D 79 (2009) 094503 [arXiv:0902.4267]. 

[7] H. Neuberger, Exactly massless quarks on the lattice, Phys. Lett. B 417 (1998) 141 
[arXiv:hep-lat/9707022]. 

[8] P. H. Ginsparg and K. G. Wilson, A Remnant of Chiral Symmetry on the Lattice, Phys. Rev. D 25 
(1982) 2649. 

[9] H. Fukaya et al. [JLQCD and TWQCD Collaborations], Two-flavor lattice QCD simulation in the 
epsilon-regime with exact chiral symmetry, Phys. Rev. Lett. 98 (2007) 172001 
[arXiv:hep-lat/0702003]; H. Fukaya et al. [JLQCD collaboration], Determination of the chiral 
condensate from 2+1-flavor lattice QCD, Phys. Rev. Lett. 104 (2010) 122002 [Erratum-ibid. 105 
(2010) 159901] [arXiv:09 11.5555]. 

[10] T. Banks and A. Casher, Chiral Symmetry Breaking in Confining Theories, Nucl. Phys. B 169 (1980) 
103. 

[11] J. J. M. Verbaarschot, The Spectrum of the QCD Dirac operator and chiral random matrix theory: 
The Threefold way, Phys. Rev. Lett. 72 (1994) 2531 [arXiv:hep-th/9401059]. 

[12] P. H. Damgaard and S. M. Nishigaki, Distribution of the k-th smallest Dirac operator eigenvalue, 
Phys. Rev. D 63 (2001) 045012 [arXiv:hep-th/00061 1 1]. 

[13] H. Fukaya et al. [JLQCD Collaboration], Lattice gauge action suppressing near-zero modes ofH(W), 
Phys. Rev. D 74 (2006) 094505 [arXiv:hep-lat/0607020]. 



7 



